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Abstract There have been some previous studies about when and through which route the alpine 
organisms ranging in the Japanese islands migrated from the continent, and how their distribution 
range has changed within the islands while adapting to the climate changes in the Quaternary. 
There are studies about alpine plants based on recent molecular information. A scenario has been 
proposed that there are two lineages. One is the lineage that had relict distribution during inter- 
glacials. The other is the lineage that migrated into the Japanese islands during later glacial peri- 
ods. On the other hand, animals have mobility and are able to move around following environmen- 
tal changes such as climate changes. Therefore, there is a possibility of having a different scenario 
from that of plants. We chose Erebia niphonica, and studied the distribution change of alpine or- 
ganisms in the Japanese islands during the Quaternary. We studied the haplotypes of base se- 
quences of 942 bp in total using №5 and СОГ. Samples are 427 individuals from the Japanese is- 
lands, Mongolia and Russia. We inferred past historical events by analyzing the haplotypes using 
the Nested Clade Phylogeographical Analysis (NCPA). As a result it has been suggested that the 
Erebia niphonica group went through isolation in more than one refugium in Hokkaido and more 
than one in Honshu during warm interglacials, and that they fragmented into different lineages 
within the Japanese islands. Also we specified the secondary contact points which support the 
NCPA which inferred that range fragmentation and dispersal were repeated. 


Key words Phylogeography, biogeography, Erebia niphonica, nested clade phylogeographical 
analysis, haplotype, ND5, COI, Quaternary. 


1. Introduction 


1.1 Phylogeographic history of the Japanese islands 


There were periods when the Japanese islands [the Japanese Archipelago] were connected 
with the Asian continent during the glacial periods and served as refugia for animals and 
plants of continental origin. During the glacial age, many of the arctic and subarctic ani- 
mals and plants migrated into the Japanese islands. During the interglacial age, distribu- 
tions of those animals and plants diminished or disappeared (Yasuda and Miyoshi (Eds), 
1998; Machida et al., 2003). There have been many studies of fauna based on fossils 
(Kawamura et al., 1989; Kawamura, 1998), and of flora based on pollen and wood fossils 
(Sohma and Tsuji, 1988; Suzuki and Takeuti, 1989; Morita, 2000). 


The major islands that constitute the Japanese islands are Sakhalin, Hokkaido, Honshu, 
Shikoku, and Kyushu. They stretch in a row from north to south (Fig. 1). There is the 
Tartar Strait (=the Mamiya Strait) between the Amur region and Sakhalin, and the Soya 
Strait between Sakhalin and Hokkaido, and the Tsugaru Strait between Hokkaido and 
Honshu, and the Tsushima Strait between the Korean Peninsula and Honshu-Kyushu. Of 
these, the Tartar Strait and the Soya Strait are shallow, and were formed later than 11,000 
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Fig. 1. The Japanese islands. Shows major islands and straits. 


years ago. On the other hand, the Tsugaru Strait and the Tsushima Strait are deep. Honshu 
separated from the Korean Peninsula and Hokkaido and became an island during the inter- 
glacial age after the Riss glaciation about 150,000 years ago. It is believed that Honshu has 
not been connected with the continent since then (Ohshima, 1990, 2000; Matsui et al., 
1998). Sakhalin and Hokkaido had been connected with the continent almost continuously 
during the period of 65,000 years from about 75,000 years ago up to about 10,000 years ago 
(Ono, 1990). They had been parts of the continent for longer periods of time than the peri- 
ods when they were isolated islands. Thus, during the glacial period the geographic condi- 
tion allowed organisms to migrate southward from the continent to Sakhalin and Hokkaido. 
On the other hand, alpine organisms that now have isolated distribution in the alpine belt of 
Honshu had expanded their distribution prior to about 150,000 years ago from the conti- 
nent, and they are the species that were affected by climate changes since then during the 
interglacial period and the last glacial period. 


The high altitude areas of central Honshu and Hokkaido in the Japanese islands have had at 
least two periods when glaciers developed, and there is mountainous permafrost in some 
limited areas in the high mountainous areas in Hokkaido and central Honshu even now 
(Iwata, 2003). Reconstruction studies of ancient flora by pollen analyses (Sakaguchi, 1989; 
Igarashi ег al., 1993; Yamada, 1998; Hoshino, 1998) suggest that the low altitude areas in 
Hokkaido were the southern limit of taiga during the last glacial period, and that there were 
forests richer with Larix gmelinii than those in the present northern Sakhalin. In the eastern 
and northern parts of Hokkaido, the distribution of taiga with Larix gmelinii as its main 
plant was divided in such a way that each of coniferous forests, grasslands, wet vegetation 
and alpine tundra formed segregated niches according to altitude and/or the geographical 
features. Thus the areas functioned as refugia for many organisms in the Far East (Ono, 
1990). 
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In Honshu there is a large rift zone called Fossa Magna that runs across the island from 
NNW to SSE at around 138 degrees of east longitude. The effect of the rift zone on the dis- 
tribution of organisms has been noted biogeographically (Su et al., 1998). However, there 
has been no study that suggests relevance of the rift zone to phylogeography of alpine or- 
ganisms. The time of its tectonic formation (it was under the sea up until several millions 
years ago) was much earlier than the time when the alpine organisms came to the Japanese 
islands. 


Regarding the understanding of the immigration routes of alpine butterflies into the 
Japanese islands and varying patterns of distribution on the islands, the limitations of re- 
gional biogeographic approaches or of studies based on the degree of specialization of ex- 
ternal characters have been pointed out (Hiura, 1971, 1977). There were, however, no alter- 
native means of approach. Now, using the new viewpoint of molecular phylogenetics, there 
have been studies on Pedicularis chamissonis (Fujii et al., 1997) and on Primula cuneifolia 
(Fujii et al., 1999). There are also phylogeographic studies on the whole Eurasian popula- 
tions of Japanese Rock Ptarmigan Lagopus mutus japonicus (Baba et al., 2001) and on the 
Small Tortoiseshell butterfly Aglais urticae (Vandewoestijne et al., 2004). 


The process of range expansion and contraction of organisms due to climate changes in the 
Quaternary is well known as an east-west migration in western Eurasia (Hewitt, 1996). But 
as the geographic characteristics of the north and south extension of the Japanese islands 
may suggest, the process in the Japanese islands is composed of north-south migration. 
Each species of alpine plants has two lineages in Honshu, northern and southern, and the 
boundary is in the northeast of the island (Fujii et al., 1997, 1999). As the primary cause of 
such a distribution character, a scenario has been suggested whereby the original population 
remained in the mountainous areas of central Honshu as refugia during the interglacial peri- 
od and later became the southern lineage, and a new lineage migrated into the Japanese is- 
lands during the glacial period after the interglacial period and formed the northern lineage 
(single refugium model). However, it could be considered to be possible that more than one 
refugium existed in the Japanese islands because of its complex mountain topography and 
the micrometeorological environments caused by it. That is to say, there could be a new 
scenario. As a result of isolation by taking refuge in more than one refugium during the 
warm periods, genetically differentiated populations have survived till now in more than one 
lineage by expanding and contracting their distribution areas adapting to later glacial cycles 
(multiple refugia model). In order to study such a possibility, it is necessary that the subject 
organisms of the study must have the following characteristics (Schmitt et al., 2006; 
DeChaine and Martin, 2004). (1) Adapting to environmental changes, they should be able 
to move to suitable surroundings. (2) Each individual does not usually move around and 
tends to stay in the maternal population, and they inherit phylogeographical information as 
a population. (3) They are able to maintain continuous populations even when the areas 
with inhabitable environment are geographically small. (4) They are currently distributed in 
wide areas even though sporadic. 


Body sizes of insects are relatively small, and they are capable of reproducing themselves 
as stable populations in relatively small habitats. Butterflies have more advantageous capa- 
bilities than plants for maintaining populations under the influence of glacial cycles, such as 
migrating to more suitable surroundings on their own. Therefore, it is highly possible that 
there is a historical scenario for alpine butterflies different from that of alpine plants. 


1.2 Systematics and ecology of the subject species 


We have chosen Erebia niphonica as the subject species of this study, and we will verify 
that this species satisfies the aforementioned requirements (1 to 4). The species is found in 
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Fig. 2. Erebia niphonica distribution ranges and sample sites; (a) Hokkaido, (b) Honshu. Pale green 
dot denotes the distribution ranges, red dot denotes the sample site. Indication unit is 7.5' 
longitude Х 5.0' latitude. Pie charts show the proportions of haplotypes within each region. 
White is private haplotypes restricted in one region. For abbreviations refer to Table 1. 


Sakhalin, the southern Kuril Islands (Kunashir Island), Hokkaido, and the alpine areas of 
Honshu, and can be classified into several subspecies (Fujioka, 1981; Inomata, 1986). In 
northeastern Eurasia the closely related Erebia neriene is distributed. There is much con- 
troversy about the taxonomic relationship of niphonica and neriene. Some taxonomists re- 
gard them as geographical variations of the same species (Fujioka, 1981; Asahi et al., 
1999), some as different species (Warren, 1936; Inomata, 1986; Gorbunov, 2001), and some 
regard only the Honshu one as an independent species (Kogure, 1979). In the monograph 
on the genus Erebia (Warren, 1936), these taxa belong to the aethiops group together with 
Е. aethiops which ranges from Europe to near Lake Baikal in southern Siberia, and Е. al- 
cmena which is found in the alpine area of central south China. The monophyly of this 
group has been verified by molecular phylogeny (Sekiguchi er al., 2002). There was a view 
that Erebia vidleri, which is found in northwestern North America and Erebia niphonica are 
closely related species (Layberry et al., 1998), but they are quite far apart based on recent 
molecular phylogeny (Nakatani et al., in preparation). 


Erebia niphonica inhabits grasslands along streams in subalpine coniferous forests above 
1,000 m above sea level in northeastern Honshu and above 1,500 m in central Honshu, 
grasslands in sparse forests of Betula ermanii near the timberline, and alpine herbaceous 
and scree grass land above the timberline (3,000-3,050 m maximum) (Nakatani and 
Kitagawa, 2001). In Hokkaido it ranges widely from near the seashore to the alpine herba- 
ceous zone around 1,600 m above sea level (Takeuchi, 2003). In Sakhalin it occurs in the 
south of the central part (N51-N52 degree), and ranges from near the seashore to grasslands 
in the forests (Asahi et al., 1999; Asahi and Kohara, 2004). 
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As for е larval food plants of Erebia niphonica, about 10 species of Gramineae апа 
Cyperaceae have been listed (Fukuda et al., 1984). In Hokkaido, larvae have also been 
found from Poa pratensis and Festuca rubra which were brought in and planted on the 
slopes of woodland trails (Kawata and Kitahara, 1990). This species is a generalist for lar- 
val food plants and appears once a year. 


The habitats in Honshu consist of relatively small fields isolated by forests, shrubs, or bam- 
boo grass lands,, and population size is small (Nakatani and Kitagawa, 2000). There has 
been no detailed study on the home range size. We have observed that the insects remained 
still in the thickets under strong wind conditions, and that they had a tendency to land 
quickly and remain still when forced to fly. It seems that usually they seldom move from 
the natal site. However, a case has been reported in which a new population started to breed 
along a new trail which had been opened in a subalpine coniferous forest (Kawata and 
Kitahara, 1990; Nakatani and Hosoya, 2003). Therefore, females seem to have some mi- 
grating habit. In other words, that the process repeats itself whereby a new colony comes 
into existence at the periphery of distribution ranges, and then becomes extinct as grassland 
shifts to shrub or bamboo grassland. 
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Fig. 3. Statistical parsimony network. Circles denote haplotypes, size of the circle reflects the 
wideness of distribution ranges in 3 levels, small black circles denote missing or theoretical 
haplotypes, solid lines connecting each pair of haplotypes represent a single mutational step, 
regardless of their length. Pie charts show the proportions of regions within each haplotype. 
Uncolored haplotype was detected in one region only. The haplotypes from the continent 
are from Khrebet Khekhtsir, Amur region, Russia (N48°15’21”, E135°03’57”) and Tsenkher 
Jiguur SPA, Hangai Mts, Mongolia (477197167, Е101%39725%. Genetic distances (%) are 
calculated in Kimura’s 2 parameter method based on СОГ (510 bp). 


We can see that the ecological characteristics of Erebia niphonica well satisfy the afore- 
mentioned necessary conditions for the study of phylogeographic history. Since the goal of 
this study is to know how the glacial cycles of the Pleistocene affected the genetic diversity 
of alpine butterflies, it is desirable to collect samples at many sites covering a wide range of 
geographical distribution (Vila et al., 2005). 


2. Materials and methods 


2.1 Specimens and collection of molecular information 


We have collected 423 individuals in 289 populations in total, including 147 individuals in 
119 populations in Hokkaido and 276 individuals in 170 populations in Honshu. Fig. 2 
shows the distribution and sampling sites of Erebia niphonica in the Japanese islands with 
7.5 longitude X 5.0” latitude as a unit, and the sampling sites cover approximately 70% of 
the distribution range. Samples were collected between years 2003 and 2005. Although 
many samples including males and females are required for reliable statistical analyses, 
from the standpoint of nature conservation, we have collected only males as samples. In 
samples from some areas, specimens were included that had been collected, dried and pre- 
served since 1996. We also analyzed 2 individuals of Erebia neriene from the Khangai 
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Бір. 4. Haplotype network and associated nested design. АП haplotypes are nested іп 4th level 
clade in both islands. 
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Table 2. AMOVA results for test of genetic subdivision between populations. 








(a) Hokkaido 

Source of variation ое % of variation 
components 

Among mountain regions 0.746* 47.48 

Among population within mountain regions 0.133 8.45 

Within population 0.693* 44.07 

(b) Honshu 

Source Of variation оше % of variation 
components 

Among mountain regions 1.492* 77.08 

Among population within mountain regions 0.165* 8.51 

Within population 0.279% 14,41 


ж Р<0.001. 


Mountains in Mongolia and 2 from the Amur region in Russia for reference. Each sample 
was caught with a net, wrapped in paraffin paper, and brought back within two days. The 
body of the sample without wings was preserved in 99% ethanol and kept at 4°C until DNA 
was extracted. 


For genetic markers, we used а part of the NADH dehydrogenase subunit 5 (ND5) and cy- 
tochrome oxidase subunit І (СОГ genes in mitochondrial DNA. DNA was extracted from 
the body and legs by DNeasy Tissue Extraction Kit (Qiagen). PCR was conducted in a 10 
ul reaction system by adding DNA (ca 40 ng), as recommended by Applied Biosystems re- 
action constitution. Thermal cycler profile was one cycle at 95°C for 10 min; 25 cycles at 
95°C for 30 sec, 45°C for 30 sec, 72°C for 1 min 30 sec, one extension cycle at 72°C for 8 
min 30 sec. PCR products were processed after Exo-SAP treatment. The nucleotide se- 
quences of the amplified DNA fragments were determined by direct sequencing method 
with Big Dye Terminator v1.1 Cycle Sequence Kit. As primers for PCR and direct se- 
quencing, we used УІ (5-ССТ СТТ TCT ССТ TTA ОТТ СА-32, C2 (5-ATC YTT WGA 
АТА AAA YCC AGC-3’) (Yagi et al., 1999) for ҮР 5, and mtD6 (C1-J-1718; 5-GGA GGA 
ТТТ GGA AAT TGA TTA СТТ CC-3’), Nancy (C1-N-2191; 5-CCC GGT AAA ATT ААА 
АТА ТАА АСТ ТС-3? (Simon еі al., 1994) for COI. We aligned the nucleotide sequence 
data with Clustal X (Thompson ег al., 1997), and conducted the fragment comparison of 
432 bp for ND5, 510 bp for COI, and 942 bp in total. 


2.2 Analysis of genetic characteristics of populations 


In order to test isolation by distance, we divided the Hokkaido population into 12 groups 
and the Honshu population into 17 groups (Table 1). We then computed the genetic vari- 
ance at each hierarchical level within a population, between populations, and between 
groups based on АМОУА (Excoffier ег al., 1992) by software Arlequin 2.0 (Schneider et al., 
2000). The number of permutations was set at 10,000 times when computing. Genetic di- 
versity of a population is considered to increase by repeating the process of range fragmen- 
tation into local groups and re-merging under climate changes (DeChaine and Martin, 
2004). We calculated standard diversity indices by Arlequin 2.0 as the indicator of genetic 
diversity of the population. 


2.3 Nested clade phylogeographical analysis 


Nested clade phylogeographical analysis (NCPA) (Templeton et al., 1995; Templeton, 
1998, 2004) can infer past events that affected the population constituents of a species such 
as past fragmentation, restricted gene flow and range expansion by combining haplotype 
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Table 3. Inferences of historical processes. 














(a) Hokkaido 

Clade Steps in inference key Inferred event х P-value 
1-5 1-2-3-4-9-Мо Allopatric fragmentaion 30.08 0.116 
1-8 1-2-3-4-9-Мо Allopatric fragmentaion 48.39 0.096 
2-2 1-2-3-4-Мо Restricted gene flow 48.49 0.047 
2-3 1-19-20-No Inadequate geographical sampling 1.00 0.021 
2-4 1-19-Мо Allopatric fragmentaion 19.00 0.000 
3-1 1-19-No Allopatric fragmentaion 45.00 0.000 
3-2 1-2-11-12-Мо Contiguous range expansion 51.98 0.000 
4-1 1-19-No Allopatric fragmentaion 127.68 0.000 
(b) Honshu 

Clade Steps in inference key Inferred event х P-value 
1-14 1-19-20-2 Inconclusive outcome* 11.00 0.002 
1-21 1-2-3-4-9-Мо Allopatric fragmentaion 117.81 0.000 
2-5 1-19-No Allopatric fragmentaion 29.00 0.000 
2-6 1-2-11-12-13-Үев Past fragmentation followed by range expansion 131.95 0.000 
2-7 1-19-20-2-11-12-13-Үев Past fragmentation followed by range expansion 1.00 0.000 
2-8 1-2-3-5-6-13-Үеѕ Past fragmentation followed by range expansion 44.81 0.000 
2-9 1-2-3-4-Мо Restricted gene flow with isolation by distance 39.44 0.000 
3-3 1-2-11-12-Мо Contiguous range expansion 82.92 0.000 
3-4 1-2-3-5-6-7-Үеѕ Dispersal but with some long distance dispersal 297.24 0.000 
4-2 1-19-Мо Allopatric fragmentaion 257.95 0.000 


x?-test of geographical association of clades and inferences of biological causes for such association. 
The inferences were obtained following the most updated key (Templeton, 2004). 
*Tip/Interior status cannot be determined. 


phylogeny and geographical information of its distribution. We inferred topology among 
haplotypes by using statistical parsimony methods (Templeton е? al., 1992) integrated into 
the TCS program (Clement et al., 2000). The network is connected with 95% reliability 
that multiple transfers do not occur at any locus. The relationships between the inferred 
phylogeny and geographic position can be evaluated by comparison of two distances of the 
within clade distance (Dc) and the nested clade distance (Dn), and the interior-tip contrasts. 
The significance of the distances defined by NCPA is evaluated using 0=0.05. These evalu- 
ations were calculated by the GeoDis Program (Posada et al., 2000). When the evaluation 
is significantly large (L) or significantly small (S), they indicate that haplotype/clade distri- 
bution is significantly widespread or significantly restricted. Significant I-T contrast for 
both Dc and Dn values denotes a displaced (L) or restricted (S) distribution of lower versus 
higher haplotypes/clades. 


3. Results 


Distribution frequency of haplotypes and sample sites is shown in Table 1. As a result, 22 
haplotypes from Hokkaido and 25 haplotypes from Honshu were detected. The number of 
substitutions of NDS is 28 (6.5%) and that of COI is 31 (6.1%). The substitution rates of 
both genetic areas are quite similar. Among 59 variable loci, those at the first codon were 
12 (20.3%), the second codon 4 (6.8%), and the third codon 43 (72.9%). Nucleotide fre- 
quencies were strongly AT biased as in many other insects (frequencies: A=30.2%, 
C=13.7%, G=13.4%, T=42.7%). There were 10 nonsynonymous substitutions, and 2 of 
them appeared in samples both from Hokkaido and from Honshu populations. 
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1st-Step Clades 2nd-Step Clades 3rd-Step Clades 4th-(Total) Cladogram 
No Dc Dn No Dc Dn No Dc Dn No Dc Dn 
1-3 
(НА001) 51.45 50.51 НС002 0.00 11.74 
НА002 0.00 37.82 НС001 0.00 11.74 


ЕТ 51.45 12.69 (НС010) 1043 11.74 
LT 6.64 0.00 












(HA023) 0.00 205.09 
(HA021) 0.00S 205.09 
(1-3) 48.205 91.535 
1-5 177.251 180.241. 
(НА003) 0.00 194.57 
НА004 0.00 194.57 
FT 131.85L 3731L 


2-1 11.745 144.46 
(2-2) 159.06 171.80 
FT 147.33L 27.34 


со 
ры 






1-5 
(НА000) 159.69 193.55 
НВ001 0.0058 191.06 
НА022 0.00 172.31 
НА020 0.00 172.31 
НА010 0.00 65.45 
FT 29.67 



























4-1 
3-1 16401L 188.30L 
3-2 93.98S 104.67S 





ЕТ not exist 
84.78 
HD001 0.00 85.92 2-3 
HD002 0.00 85.92 (1-8) 81215 81225 
HD003 21.465 50.73 (НЕ000) 0.00 122.74 
ІТ 67.40L 27.66 I-T not exist 3-2 
(2-3) 82.695 89.88 

2-4 59.345 10492 
1-12 (HF002) 26.355 72.181. I-T 23.35 -15.04 
(HF001) 53.81 59.28 1-12 48.63 54.60 
HF003 0.00 35.82 FT ~22.28 17.57 Ба 


НҒ004 31.75 29.30 
LT 37.93 26.72 





Fig. 5. Analytical results for NCPA with significant geographic structures. Nesting hierarchy pro- 
ceeds from the most recent Ist level clades on the left to the older clades on the right. 
Within each clades (boxed), columns indicate haplotypes, one hierarchical lower level 
clades or interior-tip contrasts, within clades distance (Dc), and nested clades distance (Dn). 
Individual haplotypes/clades in parentheses are interior. Significantly large and small statis- 
tics are denoted by L and S, respectively. (a) Hokkaido, (b) Honshu. 


3.1 Genetic difference among the continent, Hokkaido, and Honshu populations of E. 
neriene and Е. niphonica 


Topology of the relationships among haplotypes was inferred using a statistical parsimony 
method (Fig. 3). Nine missing or theoretical haplotypes were inferred. Due to the number 
of mutational steps separating the continental, Hokkaido and Honshu haplotypes, they 
could not be joined to the network with greater than 95% confidence. COI (510 bp) frag- 
ment comparison indicated that there is about 2.5% difference of substitutions between 
major haplotypes of the Hokkaido and Honshu populations. Also, between the Е. neriene 
population on the continent (collected in the Khangai Mountains, central Mongolia) and the 
E. niphonica population in Hokkaido or that in Honshu, there is about 2.8% and 1.8% dif- 
ference, respectively. Within the continent, the difference between those taken in the 
Khangai Mountains and those in the Amur region in Russia is only 1.0%. 


3.2 Genetic variations 


As a result of analysis of variance by AMOVA (Table 2), 48% of molecular variance in 
Hokkaido and 77% in Honshu were observed between the regions. These results suggest 
that gene flows are restricted and that there is geographic isolation between the regions. 
Standard diversity indices were calculated (Table 1). Diversity is very high in the Soya area 
(So) and in the Teshio Mountains (Te), which are located in the northern part of Hokkaido. 
In the central area (Is and Yu), diversity is low. In Honshu, diversity is very high in the Kiso 
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1st-Step Clades 


No Dc Dn 
NS131 0.00 63.26 
(NS130) 50.18 49.86 
І-Т 50.18 -13.40 
(NS150) 25.605 64.57S 
(NS121) 6.985 235.341. 
ЕТ Not exist 





94.075 
168.00L 
135.28 





NS113 0.00 
NS114 0.00 






ЕТ 90405 -57575 
(СН151) 5.39 529 
CH152 0.00 1.39 
ЕТ 539 3290 


(СН160) 2.87 


CH161 0.00 
CH162 1.04 
I-T 2.17 








CH120) 48.88 
CH121 0.005 
СН122 0.00 
СН123 0.00 53.22 
CH150 6.03S 40155 
ІТ 45.63L 6.16 


51.74L 
51.25 
53.22 





( 
















(CH130) 21.81L 21.681. 
CH131 000 9.075 
ЕТ 2181 1261 









CH141 0.00 264 
(CH140) 4.19 408 
ЕТ 4.19 145 
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2nd-Step Clades 3rd-Step Clades 
N De Dn No De Dn 


1-12 50615 5960 
(1-14) 94641 7526 
15.66 


FT 44.04L 





000 3-3 
0.00S (2-5) 64495 1758] 
16445 2-6 128.895 143.085 
97405 ЕТ  -6439S 32.74 
86.431. 

(CH110) 19.145 36.005 

1-19 4695 86.401 

LT 1445 -50405 

1-20 9605 69701 

(1-21) 49445 54225 


І-Т 39.86L -15.485 













2-9 
(1-22) 
1-23 

I-T 






20.57L 20.53L 
3.985 9.425 
16.60L 11.121 










265 


4th-(Total) Cladogram 


No De Dn 

4-2 
3-3 156.42 208.94L 
3-4 69.28S 102.525 
I-T not exist 


5b 


Mountains (Ki), the Akaishi Mountains (Ak), and the Hakusan mountain range (Ha), which 
are in the southwestern part of the distribution range. Diversity in the Joshin range (Js) is 


also high. 


3.3 Nested clade phylogeographical analysis (NCPA) 


All haplotypes were grouped into one 4th level clade (Fig.4). NCPA showed a significant 
association between genetic and geographic distribution (Fig. 5). There are 5 clades in 
Hokkaido and 9 in Honshu where a significant tendency (P<0.05) of genealogy and geo- 
graphic location was suggested by NCPA (Table 3). Clade level 1 in Hokkaido with low 
significant probability (clade 1-5, 1-8) is also shown for reference. Based on the inference 
key (Templeton, 2004), the history of Е. niphonica in the Japanese islands was punctuated 
by several cycles of allopatric fragmentation followed by range expansion. 
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3.3.1 Distribution transition in Hokkaido 
(1) Early fragmentation 


After migrating into the Japanese islands, and spreading out into areas that include the pres- 
ent distribution ranges, they dispersed into two disjunct refugia during warm interglacials, 
and divided into two lineages, the western lineage (clade 3-1) and the middle eastern line- 
age (clade 3-2, which includes the eastern and the northern lineages) (allopatric fragmenta- 
tion of clade 4-1). Then within the western lineage (clade 3-1), the populations of both 
Rishiri Island and Rebun Island (clade 2-1) were fragmented from populations in the main- 
land of Hokkaido. 


(2) Origin of northern lineage 


According to NCPA, the northern lineage (clade 2-4) is inferred to have expanded by con- 
tiguous range expansion from other clade 2-3 (these were included in clade 3-2). 


(3) Recent dispersal and habitat islanding 


There are traces that the two lineages, the eastern and the western, which had developed 
into genetically different populations during disjunct distribution in refugia, later joined to- 
gether again when they expanded their distribution in the next glacial period (see Section 
3.4). As all the populations in Hokkaido have low statistical significance at clade 1, events 
in the recent period cannot be inferred by NCPA. 


3.3.2 Distribution transition in Honshu 
(1) Early fragmentation 


A panmictic population at an early glacial period divided into two separate populations of 
the southern lineage and the northern lineage as a result of isolation in the south refugia and 
in the north refugia during the following warm interglacial period (allopatric fragmentation 
of clade 4-2). 


(2) Range expansion after the early fragmentation 


The northern lineage, which is a 2nd level clade belonging to clade 3-3, is inferred to have 
relatively long distance contiguous range expansion from the Akaishi Mts (Ak) through the 
Joetsu range (Jy) to the northernmost populations in Honshu (Chokaisan (Ch), the Mahiru 
Mts (Wa), the Hayachine range (Hy)). The southern lineage has dispersed over some long 
distance, forming the Myoko range population (clade 2-7), central Mountain population 
(the Hakusan range (Ha), the Hida Mts (Hi), the Kiso Mts (Ki), Yatsugatake (Yt)) (clade 2- 
8) and the Joshin range population (Js) (clade 2-9). 


(3) The second fragmentation, habitat islanding 


The northern population (clade 2-5) is fragmented into clade 1-13 апа clade 1-14 by al- 
lopatric fragmentation. Each of the following populations is inferred as past fragmentation 
followed by range expansion: the population (clade 2-6) which ranged from a southern part 
of the northern area to the Akaishi Mts (Ak), the population (clade 2-8) which ranged wide- 
ly over the central mountainous areas (the Hakusan range (Ha), the Hida Mts (Hi), the Kiso 
Mts (Ki), Yatsugatake (Yt)), and the Myoko range (My) population. The Joshin range (Js) 
population (clade 2-9) expanded by restricted gene flow with isolation by distance. In this 
area, there is no mountain range which constitutes a continuous geographical connection. 
Instead, relatively low mountains lie scattered separately. So, it seems a reasonable infer- 
ence. 
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(4) Recent fragmentation, multiple islanding 


Only clade 1-21 is statistically significant, and it is inferred by allopatric fragmentation with 
further fragmentation into each mountain range (the Hida Mts, the Kiso Mts, Yatsugatake). 


3.4 Secondary contacts 


A few pieces of evidence of secondary contacts can be seen in the haplotype distribution. 
As an example of the inter-lineage secondary contact in Hokkaido, a mixed area of the 
western lineage and eastern lineage is observed in Nakatoikan, Horonobecho (N44°58°37”, 
E142°02’34”) in northern Hokkaido. In Honshu, a mixed area of ће NS100 haplotype in 
the northern lineage and the CH130 haplotype in the southern lineage is observed in part of 
the Joetsu range (Jy) (N36°53°38”, E138°57’01”) (Fig. 6). In this mixed distribution area in 
the Joetsu range (Jy), mixed inhabitation is observed in all 3 populations existing within 
410 m altitudinal difference and within 2 km horizontal distance. As an example of intra- 
lineage secondary contact within the northern lineage, the haplotype NS150 which mainly 
inhabits the Yakeishi range (Yk) is observed with the haplotype NS170 together inhabiting a 
part of the Asahi Mts (As) (N44°59°55”, E142°01°36%) which lie 130 km southwest from 
the Yakeishi range. Within the southern lineage, the haplotype CH120 which ranges widely 
over the middle south area of the Hida Mts (Hi-S), Yatsugatake (Yt), the Kiso Mts (Ki) and 
the Hakusan range (Ha), is found with the haplotype CH110 in the southern area of the 
Myoko range (My) (М36%477317 E138°04’10%). At the same time, the haplotype СН110 
which ranges over the Myoko range (My) is found in the northern part of the Hida Mts (Hi- 
М) (around N36°45°42”, E137°45°32”). These areas are about 75 km apart and partitioned 
off by the Fossa Magna with 200-400 m altitude. Also the haplotype CH110 of the Myoko 
range (My) is also found together with the haplotype CH130 in the western peripheries of 
the Joshin range (Js) (N36°44°38”, E138°33°34”). Erebia niphonica is not currently found 
in the intervening areas between these secondary contact points (Fig. 7). 


4. Discussion 


Inferences of historical processes made by NCPA should be viewed with caution, because 
the inference keys make qualitative assessments, and do not give a statistical approval to a 
counter hypothesis (Knowles and Maddison, 2002). Therefore, we will examine below the 
historical transition of the distribution range of Erebia niphonica with NCPA inference as 
the basis while taking its distribution and geographical characteristics into consideration. 
We will also discuss the distribution range transition based on the ecological characteristics 
of Erebia niphonica. 


4.1 Lineage fragmentation in multiple refugia 


As can be seen in the haplotype network for NCPA (Fig. 3), there are 2 mutation steps be- 
tween the western and eastern lineages in Hokkaido, and 3 mutation steps between the 
northern lineage and southern lineage in Honshu. The first fragmentation in Hokkaido and 
Honshu can be considered to have taken place roughly at the same time. It is suggested that 
the isolation period of the two lineages in Hokkaido was relatively shorter than that in 
Honshu. This is probably because that Hokkaido was located in the north and therefore 
cooler than Honshu. 


Among the three lineages of Hokkaido, the northern lineage is separated from the other lin- 
eages by 6 mutation steps (There is one haplotype НЕ000 between them. But this was 
found in one sample among many obtained at a certain sampling sites). The isolation peri- 
od of the northern lineage is inferred to be long, but their current distribution range is re- 
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Fig. 6. The northern lineage and southern lineage of Honshu have mixed distribution areas in a part 


of the Joetsu range (Ју) (N36°53°38”, E138°5701’). 


Three mixed populations were 


confirmed to be scattered in the area within 410 m altitude difference and 2 km horizontal 
distance. Pie charts show the proportions of haplotypes within each population, orange; 
CH130 (southern lineage), blue; NS100 (northern lineage), white; NS121 (northern line- 
age). Open circles denote populations without sampled. Heavy line denotes the ridge. 
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Fig. 7. The mixed ranges of the haplotypes in northern mountainous area of central Honshu are the 
traces of the secondary contacts of the widely distributed haplotypes. Each mark denotes 
haplotype, circle; CH110, star; CH120, triangle; CH130. The shaded areas show current 
distribution ranges. Arrows show the directions of range expansions in the past. Abbrevi- 
ations denote the range names corresponding to Table 1. 
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stricted to a localized area in the north, and they coexist with the eastern lineage at many 
sites in the distribution range. Therefore, there is a possibility that they migrated into 
Hokkaido in a (later) different era from that in which the other two lineages originated. 
Recent events in Hokkaido cannot be inferred by NCPA as the statistical significance is low. 
But the topology of haplotypes (Fig. 4) suggests a possibility that many individual haplo- 
types were produced in localized areas by the bottleneck effect, while the environmental 
stresses increased during the recent warm interglacial period. 


As the lineage is fragmented into two, the northern and southern ones, in Honshu, it is sug- 
gested that there were two refugia during interglacials (Avise, 2000). We have observed in 
the alpine belt of Honshu that Erebia niphonica neither inhabits dry windward grasslands 
nor expansive grasslands without surrounding forests (Nakatani and Kitagawa, 2000, 2001; 
Nakatani and Hosoya, 2003). Considering these ecological requirements of the subject 
species, the important prerequisite for refugia is not only adequate temperature but also hu- 
midity adequate to allow coniferous forests and grasslands to grow together in patched pat- 
terns. Isolation may have also been taking place in quite dry glacial periods. This situation 
is apparently similar to the case of Erebia epiphron in Europe (Schmitt et al., 2006). 


Regarding the two lineages of Honshu, the northern lineage with the haplotype NS100 is 
scattered over a range of about 280-km between the Akaishi Mts (Ak) and the hde Mts (Ii). 
On the other hand, the southern lineage with the haplotype CH120 spreads widely from the 
Hakusan range (Ha) in the west to Yatsugatake (Yt) in the east with the Hida Mts (Hi) being 
their main distribution range. These facts suggested that the Erebia niphonica populations 
not only migrated vertically but also dispersed horizontally for long distances through 
glacials and interglacials. The dispersal seems to be phalanx type (Hewitt, 1996, 2000) as it 
occurred while keeping uniform genetic characteristics. For Parnassius smintheus in the 
Rocky Mountains, vertical movement was suggested as the main movement (DeChaine and 
Martin, 2004). This seems to be caused by the difference in the environmental require- 
ments of the subject butterflies and of the available mountain altitudes. While the habitat of 
Parnassius is a dry area at a higher altitude, that of Erebia niphonica is in the grasslands in 
subalpine coniferous forests at a lower altitude. While the Rocky Mountains are made up of 
mountains higher than 4,000 m, mountains in the Japanese islands are for the most part 
about 2,500 m high. Thus, it seems that during the warm periods Erebia niphonica was un- 
able to move to a suitable environment by moving only upward, and was forced to move 
along the latitude. 


4.2 Primary cause of vicariance of the Honshu northern lineage 


The haplotype NS100 of Erebia niphonica of the Honshu northern lineage is scattered over 
а range of about 280-km from the Akaishi Mts (Ak) to the Iide Mts (Ii) (Fig. 25). The dis- 
tribution range of Erebia niphonica in Honshu is strongly connected with the distribution 
range of Abies mariesii, which is the main constituent tree species in the subalpine conifer- 
ous forests (Nakatani and Kitagawa, 2001). In the northern area including the lide Mts (П) 
the distribution range of Abies mariesii and that of Erebia niphonica are exclusive to each 
other, and in the central area including the Akaishi Mts (Ak) they are together. The trunk of 
conifers is rather weak and breaks easily under lateral pressure. Therefore, along a moun- 
tain stream, conifers fall down under the effects of flood and avalanche. A grassland is then 
formed there. In the grassland, young conifers grow gradually, and the grassland persists 
along the stream within the forest for some time, until the young conifers grow up. 
Formation and disappearance of grassland along streams in forest take place repeatedly. 
Female Erebia niphonica from nearby colonies would visit these newly formed grasslands 
within subalpine coniferous forests, and this could result in the formation of a new popula- 
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tion. When the mountain is high enough to support an alpine herbaceous zone beyond the 
timberline, where Erebia niphonica breeds stably, it can continue to repeat expansion into 
or contraction from the grasslands forming and disappearing in coniferous forests. 
However, in low mountains where conifers cover the top of the mountain, Erebia niphonica 
cannot maintain a continuous population, even if grasslands are formed along mountain 
streams. Mountains in central Honshu are so high that they go beyond the timberline, and 
Erebia niphonica is distributed where there are conifers. But most of the mountains in 
northern Honshu are not sufficiently high to support an alpine herbaceous zone, with the re- 
sult that conifers cover the top of the mountains. Therefore, Erebia niphonica cannot occur 
there. In the northern Honshu population area, Erebia niphonica are found only in grass- 
land in the mountains where growth of conifers is suppressed by the severe winter climates, 
with heavy snow or strong wind. Mountains with altitudes of 2,500 m or higher are scat- 
tered like stepping stones in the areas from the Akaishi Mts (Ak) to the Iide Mts (Ii). These 
mountains are covered with conifers to the summits at present. But in the glacial periods, 
the timberline was lower, and there were alpine herbaceous zones at the top, in which 
Erebia niphonica could have maintained populations. In the southern Akaishi Mts (Ak) 
there are several mountains which are covered by conifers to the summits, but have grass- 
lands on steep precipices, in which Erebia niphonica is found. These populations could be 
the last survivors after the grasslands at the top disappeared. 


4.3 Genetic diversity 


In Hokkaido, the standard diversity index of each regional population is within a 0.50-0.65 
range in many areas, and regional differences in the haplotype diversity are small. 
Populations in the Soya area (So) and the Teshio Mts (Te) have high diversity. This is due 
to the endemic distribution range of the northern lineage. In Honshu high diversity popula- 
tions exist in the southwest. It is suggested that the stress of environmental changes was 
higher in the southwest, that repeated expansion and contraction of distribution ranges took 
place, and that during these processes haplotypes were fixed by the bottleneck effect. Many 
of the fixed haplotypes are different from the widely spread haplotype by 1 mutation step, 
and it may have happened in the recent warm period. On the other hand, the haplotype di- 
versity is low in the populations in the northern area. The northern area consists of relative- 
ly low mountains, and it is suggested that populations might have been reduced or even 
have disappeared by strong selection pressures during the warm period, and that there is a 
possibility that the populations in the northern area were reconstructed from small numbers 
of small populations in the next glacial period. The species that were able to survive in 
refugia in northern Honshu must have been limited to those that had high adaptability to en- 
vironmental conditions such as temperature and food plants. There are 9 alpine butterflies 
in the entire island of Honshu. But in the northern Honshu area only Erebia niphonica now 
exists. The strong adaptability of this species, which was mentioned in the Introduction, 
must have functioned effectively. 


The existence of traces of secondary contacts by different lineages, found in Hokkaido and 
in Honshu, is direct evidence that the Erebia niphonica populations underwent repeated 
fragmentation and dispersal as a result of glacial cycles. It suggests that distribution ranges 
fragmented into separate areas like islands in the past warm period, in each of which a pop- 
ulation with characteristic haplotypes evolved, and that later in the glacial period these pop- 
ulations expanded and came into secondary contact with each other, following which, in the 
next warm period the distribution ranges have again become islands. The result of NCPA 
analysis is consistent with such a scenario. These traces of secondary contacts of different 
lineages strongly support the notion that cycles of fragmentation and dispersal of distribu- 
tion ranges in the Pleistocene took place in the Japanese islands even on a small scale, just 
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as proposed for Parnassius smintheus in the Rocky Mountains (DeChaine and Martin, 
2004). In this study, the sampling locations were much more numerous than those in previ- 
ous studies and this led to the discovery of the traces of the secondary contacts. It is better 
to increase sampling locations provided that equal numbers of samples are maintained. 
This should help a better understanding of phylogeography in future studies. 


Regarding many alpine organisms in the Japanese islands, it is proposed that there are two 
lineages, the northern one which originated in the north and migrated from the Amur region 
or the Kamchatka Peninsula, and the southern one which originated in the south and migrat- 
ed through the Korean Peninsula (Saigusa, 2005). In order to take forward phylogeographi- 
cal study of animals and plants in the Japanese islands in future, a field study including de- 
tailed observations of habitats in the northern mountain area of the Korean Peninsula is 
highly desirable, as well as of those in Sakhalin, the Kuril Islands, and the Amur region in 
Russia. 
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摘 要 








第 四 紀 日 本 列島 に お ける ベニ ヒ カ ゲ の 分 布 変遷 更 (中 谷 貴 姓 ・ 宇 佐美 真一 ・ 伊 藤 建夫 ) 


日 本 列島 に 分 布 す る 高山 性 生物 が . い つど の よう な プロ セス を 経て 大 陸 か ら 渡 来 し た か , また 第 四 紀 の 
気候 変動 に 適応 し て 列島 内 で どの よう な 分 布 の 変遷 を 経た か に つい て は , 古く か ら 多 く の 考 察 が な 
され て きた が 化石 の 残り に くい 生物 に つい て の 研究 は 少な い . 最近 の 分 子 情報 に 基づく 研究 で は 高 
山 植物 エゾ コザ クラ と ヨ ツ バシ オガ マ に 関す る も の が ある . いずれ の 種 も 本 州 に 南北 2 系 統 が 存在 し 
その 境界 は 東北 地方 に ある . この よう な 分 布 状態 を 形成 し た プロ セス と し て , 南 系 統 は 本 州 中 部 山岳 地 
域 を 温暖 期 の レフ ュ ジ ア と し て 残存 し . その 後 の 寒冷 期 に 新た な 系 統 が 再度 大 陸 か ら 渡 来 し て 北 系 
統 と な っ た と いう シナ リオ が 提唱 きれ て いる ( 単 一 レフ ブフ ュ ジア モデ ル ). し か し 移動 力 を 有する 動物 で 
は , 気候 変動 な ど 環境 条件 の 変動 に よく 追随 し て 移動 で きる た め , 植物 と は 違っ た シナ リオ が 存在 する 
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可能 性 が ある . 本 州 の 複雑 な 山岳 地形 や それ に 基づく 微 気 象 な どか ら , 日 本 列島 に 複数 の レフ ュ ジ ア が 
存在 し た 可能 性 が 考え られ る . すなわち 温暖 期 に 複数 の レフ ュ ジ ア に 避難 する こと で 地理 的 に 隔離 
され た 結果 遺伝 的 に 特 化し た 個体 群 が . その 後 の 水 河 サ イク ル に 反応 し て 分 布 を 拡大 ・ 縮 小 する こ 
と で 複数 の 系 統 と し て 今日 に 至っ て いる と いう シナ リオ で ある (複数 レフ ェ ジ アモ デル ) 


この よう な 可能 性 を 調べ る た め に は , 調査 対象 生物 は 次 の よう な 特性 を 備 た て いる こと が 必要 で あ 
る . (1) 環境 の 変化 に 反応 し て 好適 な 環境 へ 容易 に 移住 で きる 移動 能力 を 備え て いる こと .(②) 個々 の 個 
体 は 通常 は あま り 移 動 を 行わ ず , 生ま れ た 集団 に と ど ま っ て いる こと で . 個体 群 と し て は 系 統 地理 的 情 
報 を 継承 し て いる こと .(③3) 生息 可能 な 環境 が 地理 的 に 狭く て も 継続 し て 生息 で きる こと .(④) 離散 的 で 
は あっ て も 現在 広域 に 分 布 し て いる こと . 


ベニ ヒビ カゲ は これ ら の 条件 を 備え て お り , 第 四 志 氷河 サイ クル の 気候 変動 に よく 追従 し て 分 布 を 拡 
大 ・ 縮 小 す る と 共に , その 歴史 的 変遷 の 履歴 を 集団 の 遺伝 的 構造 に 残し て いる 可能 性 が 高い と 期待 
され る . 


ミト コン ドリ ア DNA (mtDNA) AND5 の 一 部 (432 bp) お よび CO7 の 一 部 (510 bp) の 塩基 配列 を 決定 L, 
これ ら を 結合 し た 配列 (942 bp) に 基づい て ハプロタイプ を 検出 し た . 分 析 し た サン プル は 北海 道産 
119 個体 群 147 個体 , 本 州 産 170 個体 群 276 個体 , 合計 289 個体 群 423 個体 で ある . その 他 参 考 用 と し て 
モン ゴル ハン ガイ 山脈 産 . ロ シア アム ー ル 地域 産 各 2 個体 を 用 いた . 


検出 され た ハプロタイプ は , 北海 道産 22, AINE 25, 合計 47 タイ プ で あり , 統計 的 最 節約 法 に より 系 
統 関係 を 推定 する と , 大 陸 産 , 北海 道産 , 本 州 産 の 3 つの クレ ー ド に 分 岐 し , それ ら の 集団 間 の 塩基 置換 
は 非常 に 大 きく , 一 つの ネッ トワ ー ク 樹 に 結合 し な い (Fig. 3). サン プル 収集 地点 を 北海 道 12, 本 州 17 の 
地域 に 分 け て , AMOVA に より 分 散 分 析 を 行う と 地域 問 の 分 散 比 率 が 有意 に 大 で あり , 遺伝 子 流動 が 制 
限 され 集団 間 に 地 理 的 な 隔離 が ある こと を 示唆 し て いる (Table 2). Nested Clade Phylogeographical 
Analysis (NCPA) に より 集団 に 起こ っ た 過去 の 遺伝 的 イベ ント を 解析 し た . Inference key (Templeton, 
2004) に よる と, 日 本 列島 に お ける ベニ ヒ カ ゲ の 分 布 変遷 の 歴史 は , 異 所 的 分 断 と 連続 的 な 分 布 の 拡 
張 ま た は 分 散 の 複数 回 の 繰返し で ある こと が 示唆 され る (Table 3). 


本 州 集 団 の 4-2, 3-4, 2-8, 1-21 の 各 ク レー ド が 全て 統計 的 に 有意 と 推定 され て いる の で , これ ら に つい て 
歴史 的 プロ セス を 追っ て みる と (Table 3), 地理 的 に 分 化し て いな い 単 一 の 遺伝 的 特性 を 有する 集団 
25 温暖 期 な 間 水 期 に 異 所 的 な 2 集団 に 分 断 さ れる (クレー ド 4-2). その 内 クレ ー ド 3-3 の 集団 は 本 州 東 
北部 に 主として 分 布 す る 北方 系 統 で あり , クレ ー ド 3-4 は 中 部 地方 に 分 布 す る 南方 系 統 で ある . 次 の 
氷河 期 に 北方 系 統 (クレ ー ド 3-3) 15 range expansion, 南方 系 統 (クレ ー ド 3-4) は dispersal と 推定 され 
る よう に 両 系 統 は 共に 分 布 を 拡張 し た . 南方 系 統 で は 分 散 に 続い て , ク レー ド 2-7, 2-8,2-9 に お いて 過 
去 の 分 断 と それ に 続く 分 布 の 拡張 , ま た は 地理 的 距離 に より 制約 され た 遺伝 子 流動 が 生 じ た . ク レー ド 
1-21 で は , 引き 続い て 最新 の 分 断 が 生じ , 現在 の 分 布 域 に な みな られ る よう な 分 断 分 布 が 成立 し た . 


地理 的 に 隔離 きれ た 二 っ の 集団 間 で は 遺伝 的 に 異な る ハプロタイプ が 固定 され る . その 後 の 環境 条 
件 改善 に 適応 し て 分 布 が 拡大 する と , 再び 集団 が 出会っ て ( 二 次 的 接触 ) 二 つ の ハプロタイプ の 混 生 
状態 が 生じ る 可能 性 が ある . 筆者 ちら は 二 次 的 接触 の 痕跡 と 思わ れる 広域 分 布 型 ハプロタイプ の 混 生 
地 を 発見 する こと が で きた . すなわち 次 の 二 箇 所 で ある . (1) 北海 道 西部 系 統 の 分 布 域 で ある 幌延 町 
中 間 間 に 東部 糸 統 の ハプロタイプ (HD000) が 分 布 す る . (2) 谷川 連峰 の 一 部 で 本 州 北部 系 統 と 南部 系 
統 が 分 布 す る (Fig. б). また 同系 統 内 の 広域 分 布 型 ハ プロ タイ プ の 二 次 的 接触 跡 と し て 次 の よう な 場 
所 が ある . (3) 飛騨 山脈 の 白馬 岳 以 北 に 妙高 山系 の ハプロタイプ (CHI110) が 分 布 し , 逆 に 妙高 山系 の 
一 部 に 飛騨 山脈 の ハプロタイプ (CH120) が 分 布 する . (4) 上 信 地 域 の 岩 蒼 山 に 妙高 山系 の ハプロタイプ 
(CH110) が 分 布 す る . これ ら は 温暖 期 に 飛騨 山脈 , 妙高 山系 , 上 信 山 系 に 隔離 され た 集団 内 で , それ ぞ れ 
別 の ハプロタイプ が 固定 され , その 後 の 寒冷 期 に 分 布 を 拡大 し た 結果 二 次 的 接触 を 生じ , その 後 の 温暖 
期 に 再び 隔離 され 相手 側 の 集団 内 に 一 部 の 個体 群 が 残 存 し て いる 状態 で ある と 解釈 され る (Fig. 7). 
(5) 朝日 連峰 南西 端 に 焼 石 山系 の ハプロタイプ (148150) が 分 布 す る . 


NCPA に よる 解析 と 二 次 的 接触 の 具体 的 事例 か ら , 次 の よう な 分 布 変遷 プロ セス が 示唆 きれ る . つま 
り , 温暖 な 間 氷 期 に は 北海 道 , 本 州 と も ち に 2 ケ所 の レフ ュ ジ ア が 存在 し , 長期 に わた り 隔 離さ れ た 結 
A, 北海 道東 部 ・ 西 部 , 本 州 北 部 ・ 南 部 の それ ぞ れ 2 系 統 の 広域 分 布 型 ハ プロ タイ プ が 生じ た . ROBE 
准 期 に 分 布 を 拡大 し た 結果 二 次 的 接触 を 行い , 後 水 期 に な っ て 再び 分 断 さ れ 現在 に 至っ て いる . 二 次 的 
接触 の 痕跡 の 発見 は , 氷河 サ イク ル に 適応 し て 分 布 拡大 , 縮小 を 繰返し た 直接 的 な 証拠 で ある と 解釈 き 
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れる . 


日 本 列島 に お ける 高山 性 生物 の 由来 に つい て は , 大 陸 か ら の 複数 回 に わた る 渡来 と いう シナ リオ だ 
け で は な く , 列島 内 の 複数 地域 が レフ ュ ジ ア と し て 機能 する こと で , 複数 の 系 統 に 分 岐 し た ケー ス る も 多 
いと 思わ れる . 本 研究 で は 氷河 サイ クル に 適応 し て ベニ ヒ カ ゲ 集団 が 分 布 の 分 断 ・ 拡 大 を 繰り 返し 
た 証拠 と な る 二 次 的 接触 の 跡 を 発見 し た が , これ は 現在 の 分 布 城 た 広く カバ ー す る 地域 か ら サ ンプ 
ル を 収集 し た 成果 で あっ て , 他 の 生物 で も サン プル 収集 の 密度 を 高め れ ば , 二 次 的 接触 の 跡 が 見 つか る 
可能 性 が 高い で あろ う . 
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